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ABSTRACT 

The Feline coronavirus (FCoV) can lead to Feline infectious peritonitis (FIP), 
which the precise cause is still unknown. The theory of internal mutation suggests 
that a less virulent biotype of FCoV (FECV) would lead to another more pathogenic 
biotype (FIPV) capable of causing FIP. In this work, the 7b gene was amplified 
from 51 domestic cat plasma samples by semi-nested PCR and tested through 
phylogenetic and phylogeographical approaches. The 7b gene of Brazilian isolates 
displayed high conservation, a strong correlation between the geographic origin of 
the viral isolates and their genealogy, and its evolution was possibly shaped by a 
combination of high rates of nucleotide substitution and purifying selection. 
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76 The Feline coronavirus (FCoV) is an important pathogen of domestic and 

77 wild felids, which can cause subclinical infection, mild enteritis or lead to feline 

78 infectious peritonitis (FIP), a fatal disease characterized by inflammatory lesions of 

79 serous membranes and systemic granulomatous lesions of parenchymatous organs 

80 [23], 

81 Although the precise cause of FIP pathogenesis is still unknown, several 

82 hypotheses have been suggested [20]. The most accepted hypothesis, called internal 

83 mutation theory, suggests that during the replication of FCoV in the intestinal 

84 epithelium, a mutation occurs that makes the virus more pathogenic and able to 

85 infect monocytes and macrophages and cause FIP [23, 27], This virulent mutant 

86 variant was designated Feline infectious peritonitis virus (FIPV), while a variant that 

87 leads to enteric infection has been termed Feline enteric coronavirus (FECV) 

88 [25]. The precise nature of the mutation responsible for the pathogenesis has not been 

89 identified in the FCoV genome [10]. Nevertheless, it has been deduced that the non- 

90 structural glycoprotein 7b, codified by ORF7b, plays a determinative role in FCoV 

91 virulence [31], besides having a strong phylogenetic sign for the differentiation between 

92 FECV and FIPV [3]. 

93 To better understand the molecular epidemiology of FCoV in Brazilian domestic 

94 cats, phylogenetic hypothesis and viral population dynamics were inferred from the 7b 

95 gene. A phylogenetic hypothesis and the reconstructed population history of FCoV 

96 isolates are presented in this work, providing insights into the origins of FCoV in Brazil. 

97 Furthermore, the molecular analysis of 7b gene dispenses considerations about the 

98 internal mutation theory, regarding to the virulence of the serotypes of FCoV. 

99 This study included samples from 210 domestic cats ( Felis catus) of various 
100 breeds, random selected from different local animal hospitals (Minas Gerais, Brazil) 
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101 during 2003-2010. One hundred twenty-nine animals were healthy and taken to 

102 veterinary clinics for vaccinations and/or elective surgery. Eighty-one of them showed 

103 clinical symptoms of FIP such as anorexia, weight loss, jaundice, recurrent fever, iritis, 

104 or neurological signs and abdominal or pleural effusion [17, 23]. 

105 Blood samples were obtained by venipuncture and collected in tubes with 

106 ethylenediaminetetraacetic acid (EDTA). The plasma was obtained and frozen at -80 °C. 

107 The collection procedures were performed according to the Ethical Principles in Animal 

108 Research of the School of Veterinary Medicine of the University of Viyosa (register 

109 number 34/2010). Viral sequences isolated from healthy animals or those with clinical 

110 symptoms of FIP were designated as FECV and FIPV sequences, respectively. 

111 The complete accessory protein 7b gene (766 nt) was amplified by semi-nested 

112 PCR with two rounds of amplification using two pairs of primers previously described 

113 by Lin and others [17]. The reaction products of the semi-nested PCR were purified and 

114 sequenced by Macrogen Inc., Seoul, Korea. Contigs of the nucleotide sequences were 

115 assembled using Phred [12] and Phrap (http://www.phrap.org). The complete 7b gene 

116 coding sequences were submitted to GenBank (JX239089- JX239139). 

117 The 7b gene complete sequences of 58 FCoV isolates were downloaded from 

118 GenBank (http://www.ncbi.nlm.nih.gov/Genbank). Thus, the final dataset selected 

119 contained 109 7b gene sequences, including the sequences of Brazilian isolates. 

120 Phylogenetic evidence for recombination was tested, and recombination 

121 breakpoints were predicted using different methods (P <0.01) available in RDP3 version 

122 3.44 [19], including RDP [18], GENECONV [21], MaxChi [28], and Bootscan/Recscan 

123 [20]. Only those recombination events predicted by at least three of the methods were 

124 taken as valid; the recombinant sequences were removed from the dataset in codon 

125 selection analysis. 
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126 Selective pressure on each codon of the 7b gene sequence was evaluated using the 

127 difference between non-synonymous (dN) and synonymous (dS) substitution rates per 

128 codon using the single-likelihood ancestor counting (SLAC), fixed-effects likelihood 

129 (FEL), and internal branches fixed-effects likelihood (IFEL) methods found in 

130 DataMonkey (http://www.datamonkey.org/). 

131 Phylogenetic hypotheses for the 7b gene were inferred by Bayesian inference (BI) 

132 and maximum likelihood (ML) (Figure 1) using MrBayes v3.1.2 [16] and GARLI 2.0 

133 [34], respectively. The 7b gene sequence of canine coronavirus (GenBank ID: 

134 GUI46061) was added to the dataset as an out-group taxon to root the phylogenetic 

135 trees. 

136 Sequences were aligned using MUSCLE v.3.8.31 9 [11]. Sites with gaps were 

137 excluded. To expedite the construction of phylogenetic trees, a model of nucleotide 

138 substitution was estimated using the jModelTest program [5]. The TEVBef+I+G 

139 substitution model was selected as the best DNA evolution model according to the AIC, 

140 AICc, and BIC criteria. 

141 The BI phylogenetic trees were calculated using the Bayesian Markov Chain 

142 Monte Carlo (MCMC) method, in two runs with 50,000,000 generations and a sample 

143 frequency of 1.000. At the end of each run, the average standard deviation of the split 

144 frequencies was 0.015022. The convergence of the parameters was analyzed in 

145 TRACER v 1.5.0, and the chains reached a stationary distribution after 500,000 

146 generations. Then, a total of 1% of the trees generated was burned to produce the 

147 consensus trees. 

148 The TEVBef+Tt-G substitution model was selected in the GARLI settings 

149 (ratematrix = (0 1 2 0 3 2); statefrequencies = estimate; ratehetmodel = gamma; 

150 numratecats = 4; invariantsites = estimate), and the statistical support of the ML 
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151 phylogenetic trees was calculated by 1,000 bootstrap replicates. The 50% majority rule 

152 consensus trees of all bootstrap replicates were summarized using the SumTrees of 

153 DendroPy 3.8.0 [30]. 

154 The population history of the FCoV isolates was reconstructed using a Bayesian 

155 skyline plot (BSP), which estimates changes in the effective population size over time 

156 [8]. The BSP analysis was carried out in BEAST vl.7.2 [7] according to the BSP 

157 tutorial (http://beast-mcmc.googlecode.com/files/BSP.pdf). 

158 Only FCoV sequences were selected in BSP analysis. Sequences were aligned 

159 using MUSCLE v.3.8.31 9 [11]. Alignments were manually inspected, and the sites 

160 with gaps were excluded. The TPM3uf+I+G substitution model was selected as the best 

161 DNA evolution model by jModeltest program [8], according to the AIC, AICc, and BIC 

162 criteria. 

163 To estimate 7b gene mutation rates, the years of collection of FCoV isolates were 

164 retrieved from GenBank. Three molecular clock model assumptions (strict-clock, 

165 Bayesian-relaxed exponential molecular clock, and Bayesian-relaxed lognormal 

166 molecular clock) were tested. In each test, a MCMC run (1,000,000,000 generations) 

167 was performed considering TPM3uf+I+G as the substitution model, the respective 

168 molecular clock model assumption, and BSP as a coalescent tree prior. The high 

169 number of generations was selected to reach a large effective sample size (ESS > 200). 

170 For this purpose, analyses were processed on graphics processing units (GPUs) in a 

171 computational cluster at UFV, using BEAGLE vl.O (http://code.google.com/p/beagle- 

172 lib/) with BEAST vl.7.2. 

173 For each test, the convergence of the parameters (including the estimated mutation 

174 rate) was analyzed in TRACER v 1.5.0, and the chains reached a stationary distribution 

175 after 10,000,000 generations. The marginal likelihoods obtained in each test were 
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176 compared by Bayes factor calculations [29] with 1,000 bootstrap replicates. The test 
111 with the highest Bayes factor corresponds to the best-fit clock model and a better 

178 estimation of the mutation rate. Following this, 1% of the trees generated were burned 

179 to produce a consensus time-tree (Figure 2) using TreeAnnotator vl.7.2 [7]. 

180 To test the influence of geographic structure and of the virulence of strains in the 

181 FCoV population, the phylogenetic trees were analyzed in BaTS vl.O (Bayesian Tip- 

182 Significance testing) [22]. In these tests (geographic distribution and virulence), the 

183 high credibility set of trees estimated in the BSP MCMC run were selected, and the 

184 association index (AI) [33], parsimony score (PS) [27], and maximum monophyletic 

185 clade size (MC) [22] were calculated using 10,000 replicates (Table 1). 

186 A total of 210 plasma samples from domestic cats ( F. catus) were analyzed by 

187 semi-nested PCR from the accessory protein 7b gene. Fifty-one samples were positive 

188 for the 7b gene of FCoV. In the analysis of the positive samples was found a prevalence 

189 of asymptomatic cats of 68.63%, and 31.37% of the cats had symptoms of FIP . 

190 In sequence alignments, the 7b genes of FCoV isolates presented overall identity 

191 ranging from 41.33% (excluding sites with gaps) to 50.48% (excluding sequences with 

192 gaps). 

193 Estimation of codon selection pressures in the 7b protein showed that 27.67% of 

194 codons were predicted to be negative selection sites, with a global dN/dS estimate of 

195 0.306. Purifying selection is indicated by estimation of codon selection pressures in the 

196 7b protein. 

197 The Brazilian isolates presented higher conservation of 7b gene sequence, with an 

198 overall identity of 98.87% in the sequence alignment. Only seven polymorphic sites 

199 differentiate the sequences of JX239089 (FECV), JX239090 (FECV), JX239091 

200 (FIPV), and JX239092 (FIPV) from those of the other 47 isolates (33 FECV and 14 
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201 FIPV). These polymorphisms result in the following amino acid substitutions in the 7b 

202 protein: H160P for JX239089; H48Y for JX239090; S89F, T159N, H160P, Y167D, and 

203 C168W for JX239091; and A19S for JX239092. 

204 To describe the correlation between geographic location, virulence of strains, and 

205 genealogy estimated by Bayesian analyses, summary statistics were calculated by BaTS 

206 [22] (Tables 1 and 2) that correlate the viral phenotypic characters with the shared 

207 ancestry (represented by the phylogenetic tree). This correlation was measured by 

208 computation of the association index (AI) [33], parsimony score (PS) [27], and 

209 maximum monophyletic clade size (MC) [22]. The AI and PS test the association 

210 between traits (geographic distribution and virulence) and tree topology. The MC index 

211 tests whether traits are associated with phylogeny. Stronger phylogeny-trait 

212 associations should produce larger monophyletic clades (MC) sharing the same trait 

213 [22], 

214 The 7b gene phylogenetic trees (Figures 1 and 2) suggest a geographic pattern of 

215 the distribution of FCoV viral isolates. All Brazilian isolates (sampled between 2003 

216 and 2010) were included in the same monophyletic clade with two other North 

217 American FIPV isolates (NC_002306 and X90573, sampled in 1979 and 1981, 

218 respectively) (Fig. 1). 

219 This work provides a comprehensive analysis of the molecular epidemiology of 

220 FCoV isolates circulating in Brazil through prediction of the main events of viral 

221 introduction, and it provides new insights about viral population dynamics and selection 

222 pressures that shaped the evolution of the FCoV 7b gene. 

223 In sequence alignments, previous studies have suggested a strong correlation 

224 between insertions/deletions (indels) in the 7b gene and the virulence of FCoV viral 

225 strains [15, 31, 32]. However, no correlation between indels in the 7b gene and 
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226 virulence was found in sequence alignments, as also described by Battilani and others 

227 [2], Lin and others [17] and Bank-Wolf and others [1], 

228 With higher conservation of 7b gene sequence, the Brazilian isolates presented 

229 only seven polymorphic sites. Most of these polymorphisms were predicted to be 

230 neutral sites in the estimation of selection pressure, seeming to be random and not 

231 correlated with the virulence of strains. This high conservation of the 7b gene among 

232 Brazilian FCoV isolates shows that the internal mutation theory [1, 23, 31] possibly 

233 could not be considered for this gene. According to this theory, virulent strains (FIPV) 

234 evolve from avirulent strains (FECV) by mutation during infection in cats. 

235 Homologous recombination has an important role in FCoV evolution, and there is 

236 evidence of recombinant strains of FCoV that arose from recombination between FCoV 

237 and canine coronavirus (CCoV) [14]. In the recombination analysis of 7b sequences by 

238 RDP3, only one recombination event, involving three Taiwanese FCoV isolates 

239 sampled in 2004, was detected. Isolate DQ675437 (FECV) was predicted to be a 

240 recombinant of DQ675439 (FECV) and DQ675429 (FIPV) (P = 1.636xl0" 3 ). This 

241 finding possibly suggests a low frequency of homologous recombination of the 7b gene 

242 among FCoV strains. 

243 Substitution rates also provide good insight into virus evolution, reflecting the 

244 restrictions in genetic diversity that lead to variations in adaptability and pathogenicity 

245 of the viral population [6]. Bayes factors analysis suggested that the Bayesian-relaxed 

246 exponential molecular clock was the best-fit model for the 7b gene sequences, and the 

247 estimated mean substitution rate was 5.686xl0" 4 substitutions/site/year. This estimate 

248 agrees with what has been described for other RNA viruses, whose rates generally range 

249 from 10 2 to 10 5 substitutions/site/year [9, 13, 26]. 

250 Although most FECV and FIPV strains were included in monophyletic clades 
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251 with other viral isolates that share the same geographic origin (i.e., Brazil, the United 

252 Kingdom, the United States, or Taiwan), it was not possible to define monophyletic 

253 clades that distinguish FECV and FIPV. 

254 In geographic pattern analysis (Table 1), the topology of the 7b gene phylogenetic 

255 tree was supported by significant values of AI and PS, and all countries, with the 

256 exception of the United Kingdom (probably due to the lower sample size), showed 

257 differentiated subpopulations supported by significant MC values. In virulence pattern 

258 analysis (Table 2), no significant correlations were found by calculation of AI, PS, or 

259 MC. Thus, 7b sequences of FCoV isolates are possibly phylogenetically structured 

260 according to their geographic origin irrespective of their pathotype as shown by others 

261 [1, 2, 4, 17]. These findings contradict the hypothesis of distinct virulent (FIPV) and 

262 avirulent (FECV) strains circulating in natural populations of FCoV, proposed by 

263 Brown and others [3]. According to this hypothesis, these two viral strains would be 

264 expected to be separated into monophyletic clusters in the phylogenetic tree inferred 

265 from the 7b gene. 

266 RNA viruses may present great genetic diversity variation at the population level, 

267 allowing the reconstruction of phylogeny that reflects their epidemiological history [6]. 

268 In this way, the time tree of Bayesian Skyline analysis predicted the possible events of 

269 viral introduction over time (Fig. 2). Different events of viral introduction in Taiwan, 

270 the United Kingdom, and the USA occurred between 1850 and 1950 and are highlighted 

271 in the phylogenetic time-tree (Fig. 2). These observations are consistent with the 

272 epidemiological history of FCoV. After World War II, there was a dramatic shift in the 

273 status of cats as pets. The number of pet cats greatly increased, and this is known to 

274 favor FCoV infection [24]. 

275 A possible source of FCoV introduction in Brazil is based on the inclusion of all 
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276 Brazilian isolates in the same monophyletic clade with two North American FIPV 

277 isolates (NC_002306 and X90573), witch presented identical sequences to the 47 

278 Brazilian isolates (33 FECV and 14 FIPV) (Fig. 1). According to the phylogenetic time- 

279 tree (Fig. 2), the introduction of FCoV in Brazil possibly occurred since 1975. 

280 The authors have shown that Brazilian FCoV isolates were recently introduced 

281 from the North America, and that evolution of the 7b gene was possibly shaped by a 

282 combination of high rates of nucleotide substitution (5.686x10 4 ) and purifying 

283 selection. Furthermore, the time tree of the present study suggests that FCoV 

284 introduction in Brazil occurred over the past of 40 years. Additionally, the findings of 

285 the present study suggest that both the internal mutation theory [23, 31] and the 

286 hypothesis of distinct virulent (FIPV) and avirulent (FECV) strains circulating [3] 

287 possibly cannot be taken as valid for the 7b gene. The authors have reported high 

288 conservation among sequences of Brazilian isolates and a strong correlation between the 

289 geographic origin of viral isolates and the genealogy predicted from the 7b gene. Thus, 

290 it is more plausible that FIP is clinically manifested in cats, mainly due to host and 

291 environmental factors and independent of genetic differences between FECV and FIPV. 

292 Comparative sequence analysis may eventually not be sufficient to answer the 

293 FECV/FIPV question. 
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406 

407 FIGURES LEGENDS 

408 

409 Fig.l Evolutionary relationships between Feline coronavirus (FCoV) isolates based on 

410 the 7b gene. The majority-rule consensus tree was obtained by Bayesian MCMC 

411 coalescent analysis of 109 complete 7b gene sequences. The posterior probability values 

412 (PP) (bold; expressed as percentages) calculated using the best trees found by MrBayes 

413 are shown beside each node. The second value corresponds to the bootstrap value (BV) 

414 (underlined; expressed as percentage) that defines the clusters in the maximum 

415 likelihood tree. The outgroup taxon is an isolate of canine coronavirus (GenBank ID: 

416 GU146061). 

417 

418 Fig.2 Time tree of Feline coronavirus (FCoV) isolates reconstructed from the 7b gene. 

419 The majority-rule consensus tree of the 7b gene was obtained by a coalescent Bayesian 

420 skyline analysis with an exponential molecular clock model assumption using BEAST. 

421 The colors of branches indicate the geographic origin of the FCoV isolates 



18 


422 Figl. 



FIPV X90571 O UK 

FIPV DQ675419 © TW 

FECV DQ675421 © TW 

FECV DQ675423 © TW 

FIPV DQ675431 © TW 

FECV DQ675434 © TW 

FECV DQ675432 © TW 

FECV DQ675438 © TW 

FECV DQ675439 # TW 

FECV DQ675437 © TW 

FIPV DQ675436 © TW 

FECV DQ675447 © TW 

FIPV DQ675446 © TW 

FIPV DQ675448 © TW 

FIPV GQ152141 © TW 

FECV DQ675424 © TW 

FECV DQ675425 © TW 

FECV DQ675426 © TW 

FECV DQ675427 • TW 

FECV DQ675435 © TW 

FIPV DQ675415 © TW 

FECV DQ675433 © TW 

FECV DQ675441 © TW 

FECV FJ917522 © US 

FECV FJ943763 O US 

FIPV FJ917523 © US 

FIPV FJ917524 © US 

FIPV FJ917532 O US 

FIPV FJ917534 © US 

FIPV FJ943764 © US 

FIPV FJ943766 Q US 

FECV FJ943761 Q US 

FIPV DQ675429 © TW 

FECV DQ675442 © TW 

FECV DQ675443 © TW 

FECV DQ675444 © TW 

FECV FJ943772 Q US 

FECV JX239089 © BR 

FIPV JX239091 © BR 

FECV JX239090 © BR 

FECV JX239093 © BR 

FECV JX239094 © BR 

FECV JX239096 © BR 

FECV JX239097 © BR 

FECV JX239099 © BR 

FECV JX239100 © BR 

FECV JX239101 © BR 

FECV JX239102 © BR 

FECV JX239107 © BR 

FECV JX239108 © BR 

FECV JX239110 © BR 

FECV JX239113 © BR 

FECV JX239114 © BR 

FECV JX239116 © BR 

FECV JX239117 © BR 

FECV JX239118 © BR 

FECV JX239119 © BR 

FECV JX239121 © BR 

FECV JX239122 © BR 

FECV JX239124 © BR 

FECV JX239125 © BR 

FECV JX239127 © BR 

FECV JX239128 © BR 

FECV JX239129 © BR 

FECV 1X239130 © BR 

FECV JX239131 © BR 

FECV JX239132 © BR 

FECV JX239134 © BR 

FECV JX239135 © BR 

FECV JX239136 © BR 

FECV JX239137 © BR 

FECV JX239138 © BR 

FECV JX239139 © BR 

FIPV JX239092 © BR 

FIPV JX239095 © BR 

FIPV JX239098 © BR 

FIPV JX239103 © BR 

FIPV JX239104 © BR 

FIPV JX239105 © BR 

FIPV JX239106 © BR 

FIPV 1X239109 © BR 

FIPV JX239111 © BR 

FIPV 1X239112 © BR 

FIPV 1X239115 © BR 

FIPV 1X239120 © BR 

FIPV 1X239123 © BR 

FIPV 1X239126 © BR 

FIPV JX239133 © BR 

FIPV NC_002306Q US 
FIPV X90573 O US 

FECV F1938051 Q US 

FIPV X90578 Q US 

FIPV FJ917525 Q US 

FIPV FJ943770 © US 

FIPV F3917528 O US 

FIPV FJ917529 @ US 

FIPV DQ848678 Q UK 

FIPV FJ917527 © US 

- FECV X66718 Q US 

FIPV DQ675428 © TW 

FIPV DQ675452 © TW 

FIPV DQ675430 © TW 

FECV DQ675451 © TW 

FIPV DQ675449 © TW 

FECV X90574 Q US 

FIPV X90577 © US 

FIPV X90575 O US 

FIPV DQ675440 © TW 

FIPV DQ675450 © TW 

Canine coronavirus 


1979 

2003 

2004 
2004 
2004 
2004 
2004 
2004 
2004 
2004 

2004 

2005 
2005 
2005 

2007 
2004 
2004 
2004 
2004 
2004 

2003 

2004 
2004 

2008 
2008 
2008 

2007 

2008 
2008 
2008 
2007 

2007 
2004 
2004 
2004 

2004 

2008 

2007 
2010 

2008 
2007 
2007 
2007 
2007 
2007 
2007 

2005 
2005 
2007 
2007 
2007 
2007 
2007 
2007 
2005 

2007 
2003 
2005 

2008 
2005 
2007 
2005 
2007 
2005 
2007 
2005 
2007 
2007 
2003 
2007 
2007 

2003 

2004 
2007 
2004 

2004 

2005 
2010 

2007 
2005 

2008 

2004 

2005 
2005 
2010 

2005 

2007 

2003 
1979 
1981 
1995 
1985 

2008 
2008 
2008 
2008 

2006 
2008 
1979 

2004 

2005 

2004 

2005 
2005 
1981 
1985 
1981 

2004 

2005 


423 


0.06 Expected changes per site 



























1350 - 


19 


424 Fig2. 
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426 


Tablel. Geographic effect on the population structure of Feline coronavirus (FCoV) isolates. 


Tested correlation 

Statistics a 

Observed 

Expected b 

P-value* 

Geographic origin 

Al 

0.2413 

(0.0125; 0.5361) 

7.4867 

(6.6256; 8.2920) 

0 

Geographic origin 

PS 

7.65 

(6.0000; 9.0000) 

48.5762 

(45.6100; 51.2300) 

0 



15.11 

2.4121 


Taiwan 

MC 

(14.0000; 

23.0000) 

(2.0200; 3.2400) 

1.00E-04 

United States 

MC 

9.25 

(9.0000; 12.0000) 

1.9445 

(1.4700; 2.3500) 

1.00E-04 

United Kingdom 

MC 

1 

(1.0000; 1.0000) 

1.005 

(1.0000; 1.0200) 

1 



46.85 

3.5105 


Brazil 

MC 

(29.0000; 

51.0000) 

(2.7300; 4.5900) 

1.00E-04 


The numbers in parentheses correspond to the 95% lower and upper bounds of the highest- 
probability density intervals. 

a AI: association index; PS: parsimony score; MC: maximum monophyletic clade; 
b Expected value on null hypothesis (random phylogeny-trait association). 

^Statistical significance of tests: P <0.01. 


Table 2. Effect of virulence on the population structure of Feline coronavirus (FCoV) isolates. 


Tested correlation 

Statistics a 

Observed 

Expected b 

P-value* 

Virulence 

AI 

4.8142 

(3.7616; 5.7645) 

5.6474 

(4.7441; 6.5434) 

0.07 

Virulence 

PS 

32.03 

(29,0000; 34.0000) 

36.0344 

(32,8600; 38.87004) 

0.02 

FECV 1 

MC 

5.81 

(5.0000; 9.0000) 

4.2245 

(3.3500; 6.0200) 

0.1 

FIPV 11 

MC 

4.46 

(4.0000; 6.0000) 

3.3531 

(2.6200; 4.3200) 

0.15 


The numbers in parentheses correspond to the 95% lower and upper bounds of the highest- 
probability density intervals. 

a AI: association index; PS: parsimony score; MC: maximum monophyletic clade; 
b Expected value on null hypothesis (random phylogeny-trait association). 

I Less-pathogenic biotype of FCoV 

II Pathogenic biotype of FCoV 

* Statistical significance of tests: P <0.01. 
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